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Abstract 

A current constrained approach is proposed to calculate negative differential conductance 
in molecular nano-junctions. A four-site junction is considered where a steady-state current 
is forced by inserting only the two central sites within the circuit. The two lateral sites 
(representing e.g. dangling molecular groups) do not actively participate in transport, but 
exchange electrons with the two main sites. These auxiliary sites allow for a variable number 
of electrons within the junction, while, as required by the current constrained approach, 
the total number of electrons in the system is kept constant. We discuss the conditions 
for negative differential conductance in terms of cooperativity, variability of the number of 
electrons in the junction, and electron correlations. 



I. INTRODUCTION 



Understanding of charge transport through low dimensional electronic system has got 
tremendous impetus in recent years because of their huge potential in advanced electronic 
devices.'^ Especially, the possibility to switch a molecule between off and on states and 
to support decreasing current with increasing voltage, the so called negative differential 
conductance (NDC) behavior, are hot topics in the field of single molecule electronics, mainly 
due to its credible applications. Since its first observation in the tunnel diode by Esaki in 
1958,E1 NDC has been subject of many experimental and theoretical studies. NDC has 
been observed in a variety of experimental systemd^El g^^^ ]-^g^g ^gg^ discussed theoretically, 
even if most of theoretical work was based on one-electron pictures.'^'^ Recent experiments 
demonstrated NDC in double quantum dots junctionJ^^Hl g^nd have rekindled interest in 
the phenomenon, occurring in the low temperature, weak-coupling limit. Several theoretical 
studies on donor-acceptor double quantum dot systems address strong rectification and NDC 
features within many-electron pictures. CSHIS] However, all available theoretical studies of 
NDC are based on voltage constrained (VC) approaches, where the electric current is driven 
through the junction by imposing a finite potential drop between two applied electrodes. 
Two semi-infinite electronic reservoirs are unavoidably introduced in VC models, posing 
some fundamental problems and leading to a picture difficult to reconcile with correlated 
electrons.l^^^ 

As an alternative to VC approach, current constrained (CC) approaches have been devel- 
oped to describe transport in both meso- and molecular junctions. CC approaches avoid any 
reference to reservoirs, and quite naturally apply to correlated electrons 1 ^^ * ^^ * ^°* ^ On physi- 
cal grounds, a current can be forced in a circuit by driving a time-dependent magnetic flux 
through the circuit. The model, originally proposed by W. Kohn to describe optical conduc- 
tivity in extended systems,!^ has been applied to describe electrical transport in mesoscopic 
systems,'^^'^ and more recently in molecular j unctions .'^^'^ Alternatively, variational tech- 
niques can be adopted, and the system can be forced in a current-carrying non-equilibrium 
steady state introducing properly defined Lagrange-multipliers in the Hamiltonian. While of- 
fering an interesting alternative to the popular and successful VC methods, the CC approach 
suffers from several limitations that hinders its widespread use.l^ Not explicitly accounting 
for contacts, the CC approach relies on phenomenological models for the relaxation dynam- 
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ics to describe dissipation. In polyatomic junctions, constraints must be explicitly enforced 
to avoid charge accumulation at atomic sites within the junctions, leading to a cumbersome 
numerical problem for large and asymmetric molecular systems.'^ Finally, in VC approaches 
electrons can be exchanged between the electrodes and the molecular junction and charge 
injection/depletion phenomena can be observed upon changing the applied gate voltage. 
These phenomena are critical to many processes, including NDC, but can not be described 
in CC approaches where the current is forced through a closed circuit and the total number 
of electrons in the junctions is strictly constant. 

Here we propose a simple and effective strategy to overcome the last general problem 
inherent to CC approach. Specifically, we demonstrate that if one or more auxiliary sites 
are attached to the main current-carrying nanojunction, they can exchange electrons with 
the junction, thus representing electron source or sink. Along these lines, we present for the 
first time a calculation of NDC within a CC-based model. Describing NDC in an alternative 
picture with respect to the widespread VC-based models deepens our current understanding 
of the basic physics underlying this strongly non-linear phenomenon, as required to control 
and optimize the performance of NDC-based devices. 

II. THE MODEL AND METHOD 

Previous work from two of the present authors based on the VC approach demonstrated 
interesting nonlinear behavior in the current- voltage characteristics of two-dot systems with 
correlated electron.'^'^ Taking clue from this work, here we explore transport behavior in 
two-dots junctions within the CC approach. The study becomes interesting, especially in 
the context of the number of electrons within the junction, which plays an important role 
in the low-bias current-voltage characteristics. 

Fig. 1(a) shows a rough sketch of the model. We consider a 4 sites Hubbard model as 
described by the Hamiltonian: 

4 

Ho = eihi- ti^i+i{al^ai+i^a + h.c.) 

i=l (T=t,i 
4 

+ J2 Uifiiinii (1) 
1=1 

where (a|) annihilates (creates) an electron with spin a on site i and hi = J2a(^l(^i- 
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the sake of simplicity in the following we set constant on-site repulsions, Ui = U, and 
hopping integrals tij = t, even if different choices have been explored. Moreover we consider 
an asymmetric junction with a weak bond: —€2 = = 2, t = 0.2. 

The eigen states of Hq are stationary states and do not sustain any current. To induce 
current through the central (2-3) sites we insert these two sites into a closed circuit through 
which we impose a current according to the CC prescription: 

Hi\) = Ho- \j, (2) 

where j = —^^23 E<t('^2,o-'^3,o- ~ ^-C-) measures the current flowing through the bond between 
sites 2 and 3 (the junction region). Here and in the following e and h are adopted as 
units for charge and momentum, respectively. The field A coupled to the current enters 
the Hamiltonian as a Lagrange multiplier, whose value is fixed by the requirement that 
a finite current / =< G{X)\j\G{X) > flows through the system. The eigenstate of H{X) 
describe electrons with a finite overall linear momentum which looks rather unphysical in 
an open circuit, as described by Hq. This paradox can be solved, according to Burke et 
al. connecting the two ends of the junction (sites 2 and 3 in our model system) through 
a long, thin and ideally conducting wire so that electrons escaping from the right side of 
the junction immediately appear on the left side as to maintain fixed number of electrons 
inside the junction. As schematically shown in Fig. 1, the junction is then inserted in a 
closed loop and this allows to assign a precise physical meaning to A, that was introduced 
formally as a Lagrange multiplier. Following the seminal work of Kohn,'3Sl in fact one can 
force a current through a closed ring threading an oscillating magnetic field of frequency 
u through the center of the ring to generate a spatially uniform oscillating electric field E 
across the junction. At the leading order in the field,'^the Hamiltonian in the presence of a 
magnetic flux can be written as in Eq. [2] with A proportional to the amplitude of the vector 
potential. Quite interestingly, a similar result has been recently reported via a general and 
exact projection operator technique.'^ 

The two central sites (sites 2 and 3) are inserted in a closed circuit where the current 
is forced (cf Fig. la) and hence take active part in transport. The lateral sites (1 and 
4) represent auxiliary sites: the current does not flow through these sites, but they are 
connected to the circuit and, exchanging electrons with the main sites (2 and 3), provide 
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FIG. 1: (a) A schematic representation of the four-site system, where current flows through the 
two central sites in the closed circuit, while the lateral sites behave as sources/sinks of electrons. 
(b)The biphenyl molecule with two dangling groups, R and R', represents an example of a molecule 
that can be roughly described by the four-site model in panel (a). The red arrows describe the 
direction of current flow. 

a source/sink for the injection/removal of electrons within the circuit. Being interested 
in NDC we set U » t: the model then describes four weakly connected quantum dots. 
Alternatively, the same model may represent a minimal model for a molecular junction 
with at least one weak (poorly conjugated) bond and with two dangling groups (R and R') 
connected to the current carrying skeleton as shown, as an example in Fig. 1(b). From a 
chemical perspective, the groups R and R', modeled by the auxiliary sites of Fig. la can 
behave as electron withdrawing/donating groups with respect to the main sites. 

If the auxiliary sites (1 and 4) are disconnected from the main circuit (ti2 = ^34 = 0) the 
number of electrons participating in transport is strictly constant: having inserted the two 
central sites in an (ideal) closed circuit implies that electrons escaping from site 3 immedi- 
ately enter on 2 (or viceversa). As discussed above, the obvious constraint of constant charge 
in a closed circuit hinders the description of phenomena related to charge injection/depletion 
processes, that are instead easily captured in VC approaches where electrons are exchanged 
between the junction and the leads. However, if the auxiliary sites are connected to the 
junction the number of electrons within the circuit, n = (?t-2 + ^3)? becomes a variable: finite 
ti2 and/or t^^ allow for charge injection/depletion within the junction. 

The Hubbard Hamiltonian in Eq. (1) and its current-carrying version in Eq. (2) can 
be written on a real-space basis and the resulting matrix can be diagonalized exactly. The 
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real space basis is defined by the complete and orthonormal set of functions {(pi) that 
specify the occupation of site spin-orbit als. Both the total number of electrons, A^, and the 
z-component of the total spin, Sz, commute with the Hamiltonian, and the Hamiltonian 
matrix can be defined for a particular charge and spin sector. Here we set = 4 and work 
in the subspace relevant to the ground state, i.e., Sz=0, for a grand total of 256 basis states. 

The definition of characteristic current-voltage curves is quite subtle in CC approaches. 
As discussed above, Hamiltonian in Eq. 2 describes (at the leading order in the vector 
potential) the effect of an oscillating magnetic flux driving a current through a circuit. In 
this view, A is simply proportional to the corresponding electrical field: E oc a; A. '^^'^ However 
in the DC limit, a; — )■ 0, the driving field, and hence the driving potential vanish, leading to 
the unphysical result of a finite current at zero bias. This zero-resistance picture emerges 
because relaxation dynamics is neglected: in fact a finite potential drop is required to sustain 
current just to overcome friction in the junction. The calculation of characteristic curves in 
the CC approach then requires a model for relaxation dynamics. In particular, the current- 
carrying state, \G{X) >, is a non-equilibrium state and the system relaxes back to the ground 
state, \G{0) >. The amount of power spent to sustain the system in the non-equilibrium 
current-carrying state is set by the relaxation dynamics: the faster the relaxation the more 
power must be spent to sustain the current. The Joule law defines the relation between the 
electrical power spent on the molecule, W, and the potential drop needed to sustain the 
current: W = IV. Since / is known, characteristic I{V) curves can be obtained from 

We adopt a simple phenomenological model for relaxation dynamics that is widely used 
in spectroscopy. The model is general and properly obeys basic physical requirements, 
including detailed balance. In particular relaxation is described by introducing the density 
matrix (a) written on the basis of the eigenstates of H{0). On this basis the equilibrium 
density matrix is diagonal, with diagonal elements describing the Boltzmann populations 
of the relevant states. Relaxation dynamics of digonal elements of a are only affected by 
inelastic scattering and are written as^^ 

{o'Rjkk = y^.lkmO'mm ~ '^.ImkO'kk (3) 
m m 

where 'jkm measures the probability of the transition from m to as due to inelastic scattering 
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events. Off diagonal elements of the relaxation matrix may be written as: 

{(^Rlkm = —^kmO'km (4) 

where 

Tfcm = ilkk + 7mm) /2 + 7fcm (5) 

accounts for both inelastic scattering events, via the population inverse lifetimes 'jkk = 
J^'mlmk, and for elastic scattering events via that describes the loss of coherence due to 
pure dephasing phenomena J ^^ l ^o i 

The total work exchanged by the junction, Tr{&jiH), has 2 non- vanishing contributions: 
W = —XTr{&iij) that represents the work done on the junction to sustain the current and 
a second contribution, Wd = Tr^aj^Ho), that measures the work dissipated by the junction. 
Since the current operator, j, has vanishing diagonal elements on the real eigenstates of 
if(0), W is only affected by the relaxation of off-diagonal elements of the density matrix. 
As a result, both elastic and inelastic scattering events affect W . This is a physically relevant 
result; in fact, both phenomena concur to build up the resistance in the junction. On the 
opposite, the power dissipated by the junction, Wd = Tr^ajiHo), is only affected by inelastic 
scattering events. It represents just a fraction of the total power dissipated on the junction, 
thus suggesting that some dissipation must occur at the electrodes.'^ 

While the relaxation model is fairly crude and requires invoking electrodes to actually 
balance the energy, it has the advantage of setting on a firm basis the relation between the 
molecular resistance and dissipation phenomena. In particular we notice that the proposed 
model properly predicts the maximum conductivity of the simplest junction as e^/h as a 
result of the maximum lifetime of the electron within the junction which cannot be longer 
than the time required to the electron to cross the junction.!^ This observation drives us to 
the non-trivial role of electrodes on relaxation dynamics. In the case of very weak contacts, 
the presence of metallic surfaces close to the molecular junction is expected to open new 
channel for energy dissipation then reducing the population lifetimes.'^ On the opposite, in 
the strong contact regime the mixing of the discrete molecular eigenstates with the contin- 
uum of states of the metallic leads is responsible for increased dephasing rates.^ A very 
simple model for the relaxation matrix can then be obtained in the strong contact regime, 
where depopulation contributions to Tkm can be neglected, and in the hypothesis that the 
effects of the electrodes on the dephasing rates is similar on all states, so that = T- In 
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FIG. 2: The color map shows differential conductance (J^) as a function of potential drop (V) 



and of the energies of the auxiliary sites satisfying the relations (a) ei 
Results are shown for U=b. 



-64 and (b) ei = €4. 



this hypothesis W = XTI and the potential drop is simply proportional to A: V 
Without loss of generality we set F = 1 in the following. 



III. RESULTS AND DISCUSSION 



We consider two different choices for the energies of the auxiliary sites: in the first 
model (model a) we set ei = -64 to mimic groups with opposite electron donating/ accepting 
characteristics. In the second model (model b) we set ei = €4 mimicking two equivalent 
side-groups added to the main chain (R=R' in Fig. 1(b)). Quite interestingly, the second 
model can also describe the effects of charge injection/depletion in the junction as a result 
of an applied gate voltage. The color maps in Fig. 2 show the differential conductance ( J^) 
calculated as a function of the potential drop (V) and of the energies of the auxiliary sites, 
61 = —64 (left panel) and ei = £4 (right panel) for U = 5. NDC is not observed in either case. 
This result, confirmed at different U values, is not surprising. In fact NDC suggests a largely 
non-linear behavior of the junction, with a bistable V{I) characteristic. Such a behavior 
requires cooperativity and hence the introduction of some self-consistent interaction in the 
model. 



With this requirement in mind, we modify the model to account for possible effects of 
the voltage drop on the auxiliary sites. In fact, sites 1 and 4 do not play any active role 
in transport, but can be affected by the voltage drop at the junction. In particular, being 
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directly connected to sites 2 and 3 (device region), they can experience at least a portion of 
the electrical field required to sustain the current in the main junction. As a result, a term 
is added to the Hamiltonian affecting the auxiliary-site energies as follows: 

xV{ni — rfi), (6) 

where x measures the fraction of the potential drop that is felt by sites 1 and 4. We observe 
that only the energy of auxiliary sites must be explicitly modified by the potential: the effect 
of the potential on the junction (sites 2 and 3) is already accounted for by the Hamiltonian 
in Eq. [2j The correction term in Eq. [6] introduces cooperativity in the model: the on-site 
energies, ei and €4, do depend in fact on the current flowing through the junction that, in 
turn is affected by the on-site energies. Color maps in Fig. 3, obtained setting x = 0.5 show 
NDC regions, demonstrating the cooperative nature of the phenomenon. Fig. 3 summarizes 
the main results of this paper showing the differential conductance calculated for the two 
choices, (a) ei = —64 (left columns), and (b) ei = £4 (right columns) and U =5 and 2 (top 
and bottom panels, respectively). NDC feature is definitely more prominent in model b 
than in model a and is suppressed upon decreasing electron correlation strength, U (this 
conclusion is also supported by extensive calculations run at different U). These findings 
are in line with results obtained in the VC approach, showing that strong on-site e — e 
correlations favor NDC in double-quantum dot systems.l^ Note that, the {§^) map (and 
the I — V characteristics as well) are asymmetric because of the inherent asymmetry in the 
device region, with sites 2 and 3 having different on-site energies. 

Due to symmetry reasons, the number of electrons in the junction is constant in model 
(a), n = 2. In model (b), n is instead variable. Since the hybridization energy, t, is small 
with respect to on-site energies, n{V, ei) assumes almost integer and constant values (shown 
as italic numbers in the figures) in different regions of the V, ei plane connected by very sharp 
borders. Quite impressively, NDC features are observed for model (b) just at these borders, 
suggesting a strong relation between NDC and sharp variations of n. Upon decreasing U, the 
variation of n becomes smooth and the NDC regions become less prominent. Minor NDC 
features are also observed in model (a) with constant n, in regions where the occupancy of 
the two auxiliary sites abruptly changes from {hi — n^) =0 to ±2. 

To obtain a more detailed picture of the microscopic origin of NDC, we disentangle the 
contributions to the current from the different states. Specifically we work a real-space basis, 
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FIG. 3: Color map showing differential conductance (^) as a function of potential drop (V) and 
the energies of the auxiliary sites for ei = —64 (left panels) and ei = €4 (right panels), for two 
different U values. The blue numbers in the right panels show the approximate value of the number 
of electrons in the junction, n, in the different regions. The same number is constant (n = 2) in 
the left panels. 

where states are defined assigning electrons with a specified spin to each site in the junction. 
A generic state \p) = \0,^, i, 4t) then graphically represents the state where no electrons are 
found in the first site, one electron with spin a and one electron with spin /3 are located on 
sites 2 and 3, respectively, and site 4 is doubly occupied. In more formal language: |0, t, i, it 
) = '^2,a'^3,/3'24./3'^4,a|0) where a\^^ are the electron creation operators entering the Hamiltonian 
m eq. [l|and |0) is vacuum state- The current-carrying ground state has complex coefficients 
Cp in the real space basis, \G{X)) = J2p Cp\p) and the expectation value of the current operator 
can be written as / = J2pr ^pr, where Ipr = CpC*jpr, and jpr is an off-diagonal element of the 
current operator on the chosen basis jpr = {p\j \r) (diagonal elements of j vanish in the chosen 
basis). For the sake of clarity, we set f/ = 5 in the following discussion. Fig 4 shows the 
results relevant to model (a) with ei = —64 = 3. In this noticed before, n = 2. Two 

main contributions to / can be recognized as h = u,o, 4>lt, i> + ^\i, Ufi, t>\i, t> and 
h = I\o, n, 0, n>|o, t, 4-, n> + ^|o, n, o, n>|o, 4., t, n>- Quite interestingly, as shown in Fig. 4, 
the sum of these two contributions quantitatively matches with the total I{V) characteristic. 
The states involved in Ji and in I2 only differ by the occupation at 1st and 4th sites: the 
states involved in Ji show single occupancy of sites 1 and 4, while the states involved in I2 
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FIG. 4: Results for ei = —e^ = 3, and U = 5. The black line in the top panel shows the 
characteristic curve. The blue and green lines in the same panel show the two major contributions 
to the current as discussed in the text, the red dashed line shown the sum of the two contribution, 
and it is almost indistinguishable from the total current. In the bottom panel blue (green) lines 
show the weight of the two states responsible for the Ii {I2), as discussed in the text. The black 
dotted line is shown guide for the eyes 

have zero occupation at site 1 and double occupation at site 4. The position and height of 
the NDC peak of I2{V) exactly matches with the NDC peak of the I{V) characteristics: the 
fall of I2 and the rise of Ji at ~ 0.75 causes the NDC feature in I{V) characteristics. For 
V > 0.75, I2 starts decreasing to zero, while Ji increases to meet I2 at V"=1.04 resulting a 
valley in I{V). So, for V > 1.04, Ii dominates over I2 resulting in further increase of total 
/. 

The bottom panel of Fig. 4 shows the weight of the states that mainly contribute to 
the current, as shown in the figure label. The probability of states contributing to /i(blue 
lines) vanishes for V < 0.75: /i=0 in that region. Ii instead contributes significantly to 
total / for V > 0.75, where the contribution from I2 becomes negligible. The current 
flows only through the central bond, and, in a system with n = 2 and €2 = —£3 = 2, 
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FIG. 5: Results for ei = €4 = 2, and U = 5. Panel (a): Characteristic curve (black line) and 
total number of electrons in the junction (blue line). Panel (b) the three main contributions to the 
current in the NDC regions, as discussed in the text. The red dashed lines shows their sum to be 
compared to the total current (black line). Panel (c) shows the weight of states contributing to Ii 
(blue lines), I2 (violet lines) and I3 (green lines). More specifically the blue continuous and dashed 
line refer to states | t,i, tii) I t)it)Oii)i respectively; the green continuous and dashed lines 
refer to | |,|tit)0)) and | |, t, |t, 0), respectively; the violet continuous and dashed lines refer to 
|0,|t)tjl), and |0,t,it,l)> respectively. 

it is associated with off-diagonal elements of j mixing up states where both electrons 
reside on site 2 or are equally distributed in sites 2 and 3. NDC is observed when the 
occupancy of the auxiliary sites changes: the NDC feature occurs when the population 
of states contributing to I2, characterized by {h^ — hi) = 2, start decreasing at the ex- 
pense of the increasing population of states contributing to Ji, characterized by {h^^—hi) = 0. 

Model b is more complex, as shown in Fig. 5, collecting the results obtained for 
ei = €4 = 2 and U = 5. In this case a variable number of electrons is found in 
the junction, as shown by the blue lines in Fig. 5a. The I — V characteristic (black 
line) shows NDC features both in the negative and positive V regions, with a more 



12 



prominent feature at ~ 2. In the region of interest for NDC the current is mainly 
due to three specific contribution, as shown in Fig. 5 b. The three contributions are: 

h = n, 0, m, i, t, i) + n, o, t, i, t)j h = i\o, u, t,;>io, t, n, i) + -^lo, n, 4,t>io, i, n, t> 

and Is = J|4, t,o>|4, t, n, o)+^n, U, 4,o>|t, i, U, o>- The states involved in h have two electrons 
in the junction, n = 2, while the states contributing to I2 and Is have n = 3. As expected, 
I2 and Is contributions become sizable for large negative and positive V, respectively. How- 
ever it is clear from the figure that the NDC features observed in the characteristic I{V) 
are only due to the NDC features observed in the hiV) contribution to the total current. 
For a better understanding. Fig. 5c shows the probability of states playing a main role in 
transport. Ji drops at the NDC features because the relevant states, having two electrons in 
the junction, loose their weight and the increasing contributions of either I2 in the negative 
V region or of I3 in the positive V region, related to states with n = 3, are not able to fully 
compensate for the loss. 



IV. CONCLUSIONS 



The current constrained approach offers an alternative picture of transport in nanojunc- 
tions with respect to the widespread voltage-constrained approaches. In the CC approach 
Lagrange multipliers are used to force the system into a current carrying state, then avoiding 
any explicit reference to semi-infinite electronic reservoirs. The method then applies quite 
naturally to describe transport in systems with strongly correlated electrons. However, the 
total number of electrons taking part in transport is fixed in CC approaches, and this has 
hindered so far the application of the method to problems where charge injection/depletion 
in the junction is important. This problem could be solved introducing an effective re- 
laxation model that, accounting for the exchange of electrons between the junction and the 
electrodes, allows the system to relax towards states with different number of electrons. This 
approach however requires a fairly complex model for relaxation dynamics. Here we borrow 
from the field of optical spectroscopy a simple phenomenological model for relaxation that 
maintains constant the number of electrons, and overcome the limitation of fixed number of 
electrons connecting one or more auxiliary sites to the main junction: while the current does 
not flow through the auxiliary sites, they exchange electrons with the junction, allowing for 
a variable number of electrons taking part in transport. 
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Along these lines we attacked the problem of NDC in a very simple model with a two-site 
junction connected to two auxiliary sites. NDC is a highly non linear phenomenon: the 
bistable behavior of the V{I) curve at NDC can only be observed in largely cooperative 
systems. We introduced cooperativity imposing a dependence of the on-site energies of 
the auxiliary sites on the voltage drop: since the voltage required to support the current 
depends in turn on the energy of the auxiliary sites, the problem becomes self-consistent and 
the model supports NDC. Quite interestingly, dangling bonds (side groups) connected to the 
main current-carrying unit have recently been discussed as responsible for non-linearity in the 
characteristic curves. Particularly side groups may be responsible for quantum interference 
phenomena that drastically suppress current.'^ Our study suggests that a proper choice of 
the chemical nature of the side groups can not only reduce the current, but eventually lead 
to NDC. 

Here we just investigated an asymmetric junction with weak bonds and a grand-total of 
four electrons, in two different cases. In the first case we choose equal and opposite on-site 
energies for the two auxiliary sites leading to a constant number of electrons in the junction 
{n = 2). Minor NDC features are observed in this case, mainly related to an abrupt variation 
of the occupation of the auxiliary sites. More prominent NDC features are observed in the 
second model, where, having fixed the energies of the auxiliary sites to the same value, the 
number of electrons involved in the transport is variable. NDC is related in this case to an 
abrupt variation with the applied potential of the number of electrons in the junction. In 
all cases, large U and small t are needed to observe NDC: weak bonds and large correlations 
in fact ensure the abrupt variation in the nature of the current carrying state \G{X)), as 
required for to observe NDC. 
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